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The current emission noise of a carbon nanotube quantum dot in the Kondo regime is measured 
at frequencies f of the order or higher than the frequency associated with the Kondo effect kBTx/h, 
with Tk the Kondo temperature. The carbon nanotube is coupled via an on-chip resonant circuit 
to a quantum noise detector, a superconductor-insulator-superconductor junction. We find for 
hiy ^ ksTK a Kondo effect related singularity at a voltage bias eV ^ hu, and a strong reduction 
of this singularity for hu ~ S/cb^k, in good agreement with theory. Our experiment constitutes a 
new original tool for the investigation of the non-equilibrium dynamics of many-body phenomena 
in nanoscale devices. 
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How does a correlated quantum system react when 
probed at frequencies comparable to its intrinsic energy 
scales ? Thanks to progress in on-chip detection of 
high frequency electronic properties, exploring the non- 
equilibrium fast dynamics of correlated nanosystems is 
now accessible though delicate. In this respect, the 
Kondo effect in quantum dots is a model many-body sys- 
tem, where the spin of the dot is screened by the contacts' 
conduction electrons below the Kondo temperature Tk 
[iHS]. The Kondo effect can then be probed at a single 
spin level and in out-of-equilibrium situations. It leads 
to a strong increase of the conductance of the quantum 
dot at zero bias due to the opening of a spin degenerate 
conducting channel, the transmission of which can reach 
unity. This effect has been extensively studied by trans- 
port and noise experiments in the low frequency limit 
[4HTT] . However the noise in the high frequency limit has 
not been explored experimentally despite the fact that it 
allows to probe the system at frequencies of the order or 
smaller than kBTx/h characteristic of the Kondo effect 
[12]. In this letter we present the first high frequency 
noise measurements of a carbon nanotube quantum dot 
in the Kondo regime. We find for hv ~ ksT^ a Kondo 
effect related singularity at a voltage bias eV ~ hu^ and 
a strong reduction of it for hv ^ SksTK- These results 
are compared to recent theoretical predictions. 

The high frequency current fluctuations are mea- 
sured by coupling the carbon nanotube (CNT) to a 
quantum noise detector, a Superconductor-Insulator- 
Superconductor (SIS) junction, via a superconducting 
resonant circuit (see figure 1^). This allows us to probe 
the emission noise of the CNT at the resonance frequen- 
cies of the coupling circuit (29.5 GHz and 78 GHz) by 
measuring the photo-assisted tunneling current through 
the detector [13 . The probed sample consists of two cou- 
pled coplanar transmission lines. One line is connected 
to the ground plane via a carbon nanotube and the other 



via a superconducting tunnel junction of size 240 x 150 
nm^ (figure [l]) . Each transmission line consists of two 
sections of same length / but different widths, thus differ- 
ent characteristic impedances Zi ^ llOQ and Z2 ~ 
(figure [1^). Due to the impedance mismatch, the trans- 
mission line acts as a quarter wavelength resonator, with 
resonances at frequencies Vn = nv/Al = nz^i, with v the 
propagation velocity and n an odd integer [13]. The 
two transmission lines are close to one another to pro- 
vide a good coupling at resonance and are terminated 
by on-chip Pd resistors. The junction has a SQUID ge- 
ometry to tune its critical current with a magnetic fiux. 
The carbon nanotube (CNT) is first grown by chemical 
vapor deposition on an oxidized undoped silicon wafer 
[29]. An individual CNT is located relative to prede- 
fined markers and contacted to palladium leads using 
electron-beam lithography. The junction and the res- 
onator are then fabricated in aluminum (superconducting 
gap A = 182/ieV). A nearby side-gate allows to change 
the electrostatic state of the nanotube. The system is 
thermally anchored to the cold finger of a dilution refrig- 
erator of base temperature 20 mK and measured through 
low-pass filtered lines with a standard low frequency lock- 
in amplifier technique. 

To characterize the CNT-quantum dot, we first mea- 
sure its differential conductance dl/dVs as a function 
of dc bias voltage Vs and gate voltage Vq (figure ^) . 
For a gate voltage between 3.05 and 3.2V the CNT's 
conductance at zero bias strongly increases, a signature 
of the Kondo effect. The half width at half maximum 
(HWHM) of the Kondo ridge yields the Kondo temper- 
ature Tk = 1.4K in the center of the ridge [1 . This 
value is also consistent with the temperature dependence 
of the zero bias conductance. The Kondo temperature is 
related to the charging energy U of the CNT quantum 
dot, the coupling F to the electrodes and the position e of 
the energy level measured from the center of the Kondo 
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FIG. 1: (a) Sketch of the sample: a carbon nanotube (bottom 
electron microscope picture) is coupled to a superconductor- 
insulator-superconductor (SIS) junction (top electron micro- 
scope picture), used as a quantum detector, by a supercon- 
ducting resonant circuit. This circuit is constituted by two 
transmission lines, placed close to one another, made of alu- 
minum and terminated by on-chip Pd resistors. The SIS junc- 
tion is fabricated by shadow angle evaporation. The carbon 
nanotube is CVD grown, connected with Palladium contacts 
and side-gated, (b) Lower curve : I{V) of the detector in the 
subgap region. Upper curve: The real part of the impedance 
seen by the detector exhibits several resonances, (c) Differen- 
tial conductance dl/dVs of the carbon nanotube as a function 
of voltage bias Vs and gate voltage Vg- It exhibits a Kondo 
ridge for gate voltage between 3.05 and 3.20 V with an in- 
crease of conductance at zero bias. 



ridge, according to Bethe-Ansatz p!4| [T5]: 



Tk = VurJ^ exp 



8ur 



|4e2 - C/2 



(1) 



From U = 2.5meV, deduced from the size of the Coulomb 
diamond, and Tk = 1.4K, we obtain F = O.SlmeV. The 
asymmetry A = {Tl - ^r)/{^l + ^r) = 0.67 of the 
contacts is deduced from the zero bias conductance. 

To characterize the superconducting resonant circuit 
which couples the detector junction to the CNT, we mea- 
sure the subgap I{V) characteristic of the junction which 
depends on the impedance of its electromagnetic environ- 
ment [17 . In the case of a superconducting transmission 
line resonator [131 HI] ? resonances appear in the subgap 
region Vd < 2A/e due to the excitation of the resonator 
modes by the ac Josephson effect [18 . These resonances 
are related to the real part of the impedance Z{u) seen 
by the junction: 



I{Vd) = Re[Z{2eVD/h)] I^/2Vd , 



(2) 



with Ic = 7rA/(2ei?Ar) the critical current [18], Rn = 
28.6kQ the normal state resistance of the junction and 
A the superconducting gap of the electrodes. Equa- 
tion |2] accounts for the effect of the electromagnetic 
environment on the tunneling of Cooper pairs through 



the Josephson junction \TP . Figure [T]3 shows the I{V) 
of the junction in the subgap region for Ic maximized 
with magnetic flux. The subgap resonances thus yield 
via equation |2] Re[Z{u)] (Fig. ^p) which is peaked at 
frequencies ui = 29.5 and = 78GHz. Using the 
height and width of the resonance peaks of Re[Z{i')], 
we infer the coupling between the junction and the 
CNT [13]. We then translate a photo-assisted tunneling 
(PAT) quasi-particles current measurement into a cur- 
rent emission noise measurement for the frequencies ui 
and The ratio between the measured PAT cur- 
rent through the detector and the current emission noise 
of the CNT at a given resonance frequency is estimated 
as follows, rn is given at the resonant frequency Un by 
e^\Zt{jyn)\^Sun/{hUn)^lQP^o{VD + ^) with the tran- 
simpedance of the coupling circuit, defined as the ra- 
tio between the voltage fluctuations across the detector 
and the current fluctuations through the source, Svn the 
width of the resonance peak and Iqp^o{Vd) the I-V char- 
acteristic of the detector [13 . This value has been cali- 
brated in a previous experiment with the same design as 
the one used in the present work [13]. is then calcu- 
lated using the ratio of the calibrated sample corrected 
according to the square area under the corresponding 
peak of Re[Z{u)] (figure [T]3) , the value of the super- 
conducting gap and the tunnel resistance of the detector 
junction. 

To measure the quantum noise of the CNT, we modu- 
late its bias voltage Vs and monitor the modulated part 
of the PAT current through the detector for a given de- 
tector bias voltage Vd- Vd selects the frequency range 
of the measurement [13 . We have thus access to the 
derivative of the PAT current versus CNT bias voltage 
dIpAT / dVs at a given frequency. Using the previously 
estimated coupling coefficient, we translate this quan- 
tity into the derivative of the current noise Si at one 
of the resonance frequencies versus V^, dSi/dVs- This 
quantity is plotted in the center of the Kondo ridge, i.e. 
e = at two frequencies (Figure [2|i and b). For each fre- 
quency, the data exhibit a region close io Vs = ^ where 
dSi/dVs = 0. This corresponds to \eVs\ < hv^ where 
the system does not have enough energy to emit noise at 
a frequency v. The observation of this zero noise region 
is a strong evidence that we are indeed only measuring 
the emission noise of the CNT. For \eVs\ > hu the sys- 
tem emits noise at u. For the first resonance frequency 
ui = 29.5GHz, with hui ~ kpTx^ the measured deriva- 
tive of the noise shows a singularity for bias voltages 
close to the measured frequency. At higher bias volt- 
ages dSj/dVs is much smoother. For hvs ^ 2.1kBTK the 
previous singularity is nearly absent and dSi / dVs versus 
Vs is practically flat. 

The high frequency noise of quantum dots in the 
Kondo regime has been studied theoretically at equilib- 
rium using the numerical renormalization group (NRG) 
technique [19 . Non-equilibrium results for the finite- 
frequency noise are theoretically much more demand- 
ing. They were obtained only for peculiar values of 



3 




FIG. 2: (a) Derivative of the current noise and differential 
conductance of the CNT versus Vs in the center of the Kondo 
ridge. Left axis : measured dSi/dVs (in red) at 29.5GHz as 
a function of bias voltage Vs- The black solid line corre- 
sponds to the calculated dSi/dVs with a fitted decoherence 
rate incorporating both intrinsic and extrinsic decoherence 
(see text), while the dashed line is the theoretical prediction 
with the intrinsic decoherence rate. Right axis : Differential 
conductance dl/dVs of the CNT in units of /h. (b) Same 
data at 78GHz. Inset : measured conductance, symmetrized 
with respect to the bias voltage, in units of /h. The black 
solid line correspond the the FRG calculation of the conduc- 
tance with the same fitted decoherence rate as the one used 
for the noise measurements, (c) Schematic picture of a quan- 
tum dot in the out-of-equilibrium Kondo regime together with 
the density of states for two distinct bias voltages. When the 
quantum dot is voltage biased, the Kondo resonance splits in 
two with a splitting given by the applied bias voltage. This 
leads to an increase of the emission at frequency eVs = hu. 
The amplitude of the resonance peaks, and thus the emission 
noise at eVs = hu, can be reduced by decoherence effects 
induced by the applied bias voltage. 

parameters (strongly anisotropic exchange couplings) of 
the Kondo problem using bosonization methods [20] , 
and by using non-equilibrium real time renormalization 
group approaches [21] [22j. The latter approaches as- 
sume hv^eVs ^ ksT^^, with the Kondo temper- 
ature defined from the renormalization group. Impor- 
tantly, differs from Tk (defined experimentally as 
the HWHM of the differential conductance) by a nu- 
merical factor, which has to be determined (see below). 
Here we employ the real time functional renormalization 
group (FRG) approach developed in Ref. [22 to com- 
pute the non-equilibrium frequency-dependent noise and 
compare it to the experimental results. We perform the 
non-equilibrium calculations using the Kondo Hamilto- 
nian, given by : 

a,(3=L,R a, a' 

Here the ja(3 denote the Kondo couplings, a, P are indices 



for the left (L) and right (R) leads, a stands for the three 
Pauli matrices, and the operator i/jaa destroys an electron 
of spin a in lead a G {L, R}. We parametrize the dimen- 
sionless exchange couplings ja(3 as ja(3 = j VaV(3^ with 
the factors {vl^vr} = {cos(0/2), sin(0/2)} accounting 
for the asymmetry of the quantum dot, cos(^) = A, and 
(j) related to the T = conductance as G{T = 0) = 
{2e^/h) sin^((^). 

The Kondo Hamiltonian assumes that charge fluctu- 
ations in the CNT quantum dot are frozen. There- 
fore, the theoretical results based on ([3| can and shall 
be compared with experimental ones only for bias volt- 
ages Vs <^ U/e. As a first step, to determine the ratio 
T^^/Tk^ we computed the equilibrium conductance by 
using NRG [30 and compared it to experimental data in 
the center of the Kondo ridge (e = 0). This enabled us 
to establish Tk ~ 3.7 (see the appendix). There- 
fore, the condition hv ksT^^ for our FRG approach 
to apply is certainly met for the frequency z/3, and still 
reasonably satisfied for ui. 

Within the Kondo model, we can express the Fourier 
transform of the emission noise Sj as : 

(4) 

where 5 is a dimensionless function, which we calculate 
by solving numerically the FRG equation (see appendix) . 
Since the measurement temperature satisfies T <C T^, 
we have taken T = in the calculations and checked 
that a finite but small temperature does not affect our 
results. Note that no fitting parameter has been in- 
cluded at this level, since the asymmetry parameter A 
and ^ 0.38 K were extracted from the experimen- 
tal data. The dashed lines in figure [2^ and b show the 
calculated dSj/dVs curves for frequencies ui = 29.5GHz 
and Us = 78GHz, respectively. The computed curves 
are only shown in the bias range \Vs\ < 1 mV, where 
the Kondo Hamiltonian in equation ([3| is appropriate to 
describe the physics of the CNT quantum dot. For both 
frequencies, the theoretical curves exhibit sharp singular- 
ities at eVs = hu^ much more pronounced than the ex- 
perimental ones. This especially holds for the resonance 
frequency, us = 78GHz, where the resonance is almost 
completely absent experimentally. The singularity at the 
threshold eVs ~ hu is related to the existence of two 
Kondo resonances associated with the Fermi levels of the 
two contacts. Inelastic transitions between them lead to 
an increase of dSj/dVs for frequencies corresponding to 
the energy separation, hu ^ eVs (figure [2]3). 

To compute the dashed curves in figure [2j an in- 
trinsic spin decoherence time rs induced by the large 
bias was included and calculated self-consistently in the 
FRG approach [23 (see appendix). The decoherence of 
the Kondo effect induced by a large dc voltage bias is 
a well-known feature which has indeed been measured 
[24| [25] . and has been predicted to lead to a strong 
reduction of the Kondo resonance due to inelastic pro- 
cesses [23l [26l [27] . Since the singularity in the noise at 



4 




3.00 3.10 3.20 



Vg(V) 



FIG. 3: F{Vs) = [dSi/dVs]/[edI/dVs{Vs - hu/e)] as a 
function of the bias voltage Vs and the gate voltage Vg at 
z/i = 29.5GHz (a) and = 78GHz (b). F{Vs) is arbitrarily 
fixed to zero for e\Vs\ < hu. The gray curves on top of the 
colorplot correspond to the bias voltage dependence of F{Vs) 
at Vg = 3.12V and Vg = 3.21V. The black arrow indicates 
a value of the F(Vs) equal to 1. When the conductance is 
low, F{Vs) is close to one while it is reduced in the highly 
conducting regions. 



hu ~ eVs is associated with the transitions between the 
two Kondo resonances pinned at the Fermi levels of the 
contacts, this singularity is also affected by decoherence. 

However, as shown in figure [2j the computed intrin- 
sic decoherence time is insufficient to explain the exper- 
imentally observed suppression of the peak in dSi/dVs- 
Therefore, we incorporated a voltage-dependent spin re- 
laxation rate in our calculations, r^^(Vs'), which includes 
external decoherence. The consistency of this approach 
can be checked against the experiments: a single choice 
of T^^iVs) must simultaneously reproduce the voltage 
dependence of the differential conductance through the 
dot dI/dVs{Vs), and those of the vi = 29.5 GHz and 
1^1 = 78 GHz noise spectra, dSj /dVs{Vs). Further- 
more, should be suppressed for Vs < . We 
found that a bias-dependent decoherence rate of the form 



h/rs ~ a kBT^^aidin{/3eVs/kBT^^) (similar in shape 
to the calculated intrinsic spin relaxation rates), with 
a = 14 and P = 0.15 satisfied all criteria above. The 
continuous lines in figure [2] show the dSj/dVs curves 
computed with this form of h/rs^ and fit fairly well the 
experimental data for both resonator frequencies. As a 
final consistency check, we also computed the differen- 
tial conductance through the dot (taking into account 
the above form of r^^) and compared it to the measured 
dl/dVs curves. A very good agreement is found with- 
out any other adjustable parameter in the voltage-range 
V > O.lmF, where the FRG approach is appropriate 
(Inset of figure ^) . 

From the theoretical fits we infer that the experimen- 
tally observed noise spectra and differential conductance 
can be understood in terms of a decoherence rate, which 
is about a factor of 2 larger than the theoretically com- 
puted intrinsic rate (see appendix). One possibility for 
this discrepancy is that the experimentally observed de- 
coherence is intrinsic, and FRG - which is a perturba- 
tive approach - underestimates the spin relaxation rate 
in this regime (which is indeed almost out of the range 
of perturbation theory). Another possibility is that the 
experimental set-up leads to additional decoherence. 

The experiment also allows to draw a complete map 
of the noise in the region of the Kondo ridge. We define 
F{Vs) = [dSi/dVs{Vs)]/[edI/dVsiVs - hv/e)], i.e. the 
ratio of the derivative of the noise to the differential con- 
ductance shifted in voltage by an amount corresponding 
to the measured frequency. For both linear and non lin- 
ear systems with energy independent transmission at low 
temperature this quantity is equal to the Fano factor [28 . 
We have plotted F{Vs) for ui = 29.5GHz (figure |3[i) and 
Us = 78GHz (figure [sJd). For leVs-l < hu^ where the emis- 
sion noise is zero, FjVs) is arbitrarily fixed to zero. For 
both frequencies the noise is found to be sub-poissonian, 
with F{Vs) close to one in the poorly conducting regions 
and a strong decrease of F{Vs) along the conducting re- 
gions. This is qualitatively consistent with the reduction 
of the Fano factor for a conducting channels of transmis- 
sion close to one. This result has to be contrasted with 
back scattering noise measurements in the Kondo regime 
at low frequency and low bias voltage where the Fano 
factor was found to be higher than one |TOl [11] . 

In conclusion we have measured the high frequency 
current fluctuations of a carbon nanotube quantum dot 
in the Kondo regime by coupling it to a quantum de- 
tector via a superconducting resonant circuit. We find 
that the noise exhibits strong resonances when the volt- 
age bias is of the order of the measurement frequency in 
good agreement with theory provided that an additional 
decoherence rate is included which prevents the full for- 
mation of the out of equilibrium Kondo resonances. Our 
experiment constitutes a new original tool for the inves- 
tigation of the non-equilibrium dynamics of many-body 
phenomena in nanodevices. 
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Appendix : Relation between Tk and T^^ 

To perform the theoretical calculations, one first needs 
to find the Kondo temperature. However, the Kondo 
temperature is defined only up to a prefactor, and its 
value also depends slightly on the physical quantity 
from which it is defined. Our experimental Kondo tem- 
perature, Tx, is defined as Tk = e AVs/ks^ with 
AVs the half-width at half maximum of the measured 
G{Vs) = dl/dVs curves. In the FRG calculations, on 
the other hand, it is defined as a scale (frequency), 
^bT^^ = hu^^, where the so-called leading logarithmic 
calculations yield a divergent interaction vertex at T = 
temperature. It can, however, also be defined as the tem- 
perature, T^~^^^, at which the linear conductance drops 
to half of its T = temperature value. The ratios of 
all these Kondo temperatures are just universal numbers 
(apart from a possible but presumably small dependence 
of AVs- on the anisotropy. A). 

To determine the ratio of and ksT^^^ 

we performed numerical renormalization group calcula- 
tions [31 : We computed the full G(T) curve, extracted 
from it the width T^~^®^, and for the same parame- 
ters, we also computed = Tiuj^ /ks from the high- 
frequency tail of the so-called composite fermions' spec- 
tral function, scaling as ^ C/\ii'{(jo/uj^) at large fre- 
quencies. In this way, we obtained a ratio 



K 



-dep 



0.3. 



The ratio tJ. '^'"^ /Tk was then determined from experi- 
mental data [32 , giving 



T^T-dep 



K 



1.2 . 



The previous equations yield the ratio, /Tk ~ 0.36, 
used in our calculations and quoted in the main text. 



Appendix : Summary of the functional 
renormalization group approach 

In this work we used the functional renormalization 
group approach developed in Ref. [22] , an extension of the 
formalism of Ref. ^3^. In this approach, formulated at 
the level of non-equilibrium action, a short time cut-off a 
is introduced, and increased in course of the renormaliza- 
tion group (RG) procedure to eliminate the high-energy 
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FIG. 4: Voltage dependence of the fitted and intrinsic relax- 
ation rates, and their difference, termed as extrinsic spin re- 
laxation rate in the voltage range, Vs > 0.1 mV. The intrinsic 
rate was calculated by FRG, while the form and overall mag- 
nitude of the extrinsic rate were inferred from the differential 
conductance data. 



degrees of freedom. This procedure yields a retarded in- 
teraction (ja.^ ga,^(t — t'^d)\ whose cut-off dependence 
is described by the differential equation. 



dg(^, a) 
din a 



= g(^,«) q(^,«) g(^,«) • 



(5) 



Here we introduced the matrix notation, ga,^{uj^d) 
g(a;, a) for the Fourier transform of the retarded interac- 
tion, and the matrix q(cj, a) denotes a cut-off function. 
In our calculations we have not approximated this lat- 
ter with 6-functions, as in Ref. [22 , but used a function 
corresponding to the real time propagators of Ref. [22] , 
also incorporating the effect of an exponential decay rate 
r = l/rs. 

The voltage-dependent decay rate T{Vs) has an intrin- 
sic part, rintr(Ks')7 as well as an external contribution, 
rext(V5'). The former contribution can be identified as 
the Korringa spin relaxation rate, and can be expressed 
as 

Fintr = ^ / dujgaf3{(^) g(3a{^) 

f{uj - fla) (1 - /(^ - M/?)) , (6) 

with ga(3{^) the vertex functions in the limit a ^ oo, 
/j^a the electro-chemical potentials of the leads, and / 
the Fermi function. The rate Fintr CKsj^ as computed 
by FRG is shown in Fig. [4j Rather surprisingly, in the 
cross-over regime, eVs ~ ksT^^^^, it almost saturates, 
and only weakly depends on the bias. For a compari- 
son, figure [4] also shows the total fitted relaxation rate, 
F(Vs') = ^/ts{Vs) needed to reproduce the differential 
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conductance and noise data. It is not very far from the 
calculated intrinsic contribution, but it is above the lat- 
ter, and it apparently includes some extrinsic spin relax- 
ation, too. 

In the formalism of Ref. [22| , the current operator and 
the current vertex are also renormalized during the RG 
procedure, and also become non-local. However, the cur- 
rent vertex, V, has a more complicated structure than 
the interaction vertex, and possesses two non-trivial time 
arguments, ga^{^) ^ V'q;/3(^i, ^2)- The evolution of V 
under the RG is described by a differential equation simi- 
lar to equation ([5| (see reference [22] for the details). The 
noise spectrum, i.e., the Fourier transform of the current- 
current correlation function can then be expressed as a 



double integral of this retarded current vertex taken in 
the limit a ^ 00. The emission noise, Se{oo)^ e.g., can be 
expressed as 

5eH = ^ 5(5 + 1) 1^ Tr{V(a;+,a;_)G>(a;_) 
\{C:.,u:+)G<{u:+)} , (7) 

with cDi = a; ± ^ and S = 1/2. Here the trace refers to 
the labels = L/R^ and the bigger and lesser Green's 

functions are given as G^^^(co') = ±i 27r da/s f{^{^ — 

Ma)). 
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